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Abstract 

We use the recently introduced theory for the kinetics of systems of classical particles 
to investigate systems driven by Smoluchowski dynamics. We investigate the existence of 
ergodic-nonergodic (ENE) transitions near the liquid-glass transition. We develop a self- 
consistent perturbation theory in terms of an effective two-body potential and work to 
second order in this potential. At second order, we have an explicit relationship between the 
static structure factor and the effective potential and choose the static structure factor in the 
case of hard spheres to be given by the solution of the Percus-Yevick approximation for hard 
spheres. Then, using the analytically determined ENE equation for the ergodicity function, 
we find an ENE transition for packing fraction rj greater than a critical value rj* = 0.76 
which is physically unaccessible. The existence of a linear fluctuation-dissipation theorem 
in the problem is shown and used to great advantage. 
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I. INTRODUCTION 

We continue^ our presentation of a self-consistent approach to the kinetics of classi- 
cal systems of particles by studying the fluctuations in equilibrium of a system driven by 
Smoluchowski dynamics. We first show that the system obeys a linear fluctuation-dissipation 
theorem (FDT). This simplifies the structure of the theory significantly giving the conven- 
tional linear relation between the density-density correlation function and the conjugate 
response function. We present here a perturbation theory valid to second order in an ef- 
fective interaction potential. Because of the self-consistency we are able to show that this 
expansion is useful even for systems with hard-core interactions. 

This approach was demonstrated at first order in Ref. |l| (henceforth referred to as 
FTSPD) where the effective potential is found to be proportional to the direct correla- 
tion function. Here we extend the calculation to second order. In this paper we focus on the 
collective (or Ornstein-Zernike) second order self-energies. We show that these self-energies 
are quadratic functionals of the full density- density correlation functions, the components 
of the self-energies individually satisfy a FDT, no wave-number or frequency cut offs are 
needed 2 , and the set of Dyson equations fundamental to the theory can be replaced by a 
single-kinetic equation of the same form as that produced in memory-function theories^. 

Going further, one can show that the collective contribution, again using the FDT, gives 
a relation for the static structure which agrees with the result from a purely static calcu- 
lation giving the structure factor in terms of the direct correlation function and the direct 
correlation function self-consistently in terms of the static structure factor. Our approach is 
to assume the static structure factor is known and to solve for that effective potential which 
satisfies the second order structural equation. This effective potential is then available in 
the dynamical calculations. It could also be used in other static calculations. We focus here 
primarily on the large wavenumber regime near the structure factor maximum. 

We also show that it is only the collective part of the self-energy which enters into the 
determination of the ergodic-nonergodic (ENE) phase-diagram. The ENE transition occurs 
when the density- density correlation develops a zero-frequency 5 function. The amplitude of 
this 5 function satisfies a nonergodicity equation. This equation is very similar to that found 5 
in mode-coupling theory (MCT). For a system characterized by a structure factor obtained 
as the solution of the Percus-Yevick approximation^ for hard-spheres we find a transition 
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at packing fraction vf = 0.76. Conventional MCT gives rj = 0.51. As obtained here, rf is 
well above any physically attainable density in agreement with experiment and simulation. 
If we drop self-consistency and use the first order effective potential in determining if we 
obtain the value vf = 0.60 even though the first and second order effective potentials are 
very similar. Therefore rj* is a rather sensitive quantity. Despite our prejudices that there 
is no physical ENE transition in the single- component hard-sphere Smoluchowski dynamics, 
more work checking perturbations to our solution is required before one can claim a "proof 
of the lack of a transition. 

Experience guides us to be skeptical 7 about the stability of this transition under the 
inclusion of higher order terms in perturbation theory. Happily one can test this skepticism 
by going to higher order in the expansion. This appears to be technically feasible. 

In a companion paper we discuss the single-particle contribution to the second order self 
energies. This term in the second order self-energy makes no contribution to the structural 
statics (static structure factor) but does involve the equation of state. It also does not 
come into the nonergodicity-equation determining the ENE phase diagram. This does not 
mean that it plays no role in the slow dynamics of the system. Quite to the contrary, this 
term depends linearly on G pp and in a way which suggests the F 12 model of Goetze^ and a 
mechanism for stretching the dynamics. 

We find that the self-energy Tbb is a functional, to arbitrary order in the effective po- 
tential, of the full-density-density correlation function. Setting aside the difficult question 
of vertex renormalization, the construction of Tbb to arbitrary order appears feasible. This 
involves construction of the self-energy Tbb as a polynomial in G pp which produces a zero- 
frequency 5 function in response to the a zero-frequency 5 function developing in G pp in the 
nonergodic state. 

We are interested in the kinetics of colloidal suspensions. They are of interest because one 
can apparently carry out clean experiments in the regime where the system becomes glassy 
or crystalizes^. The system is also of interest because the colloidal particles are buffeted by a 
bath of smaller particles which has the effect of rapidly thermalizing the momentum degrees 
of freedom. This process is carefully described by Fokker-Planck (FP) dynamics. In FP 
dynamics the Newtonian equations of motion are supplemented by noise terms. It is very 
convenient to study a simpler dynamics than Fokker-Planck or Newtonian dynamics^. To a 
first approximation in colloids we can assume that the momenta thermalize quicker than the 
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position variables. Smoluchowski dynamics^ 1 assumes that the momenta are thermalized 
and particles move via a random walk or diffusion process interacting with the same two- 
body forces as in the Newtonian case. In the noninteracting limit one has iV-Brownian^ 
particles. 

While the Smoluchowski dynamics offer a complete self-consistent dynamical system with 
a static structure equivalent to that for Newtonian and FP dynamics, there has been a search 
for simpler, coarse-grained field theories applicable to the colloidal system. Considerable 
energy has been focused on the Dean-Kawasaki (SDK) model^ 1 ^. The MSR action for this 
model 15 is given by 



A nK = / d\ 



Dp{Vpf + ip 



dp / „5F\ 
-V ^ 1 



(1) 



dt V bp) 

where p is a density continuous field and p is its response conjugate. The functional derivative 
of the effective free energy, taken from density functional theory, is given by 

5F r 
— = T In p(x) + J d d yu{x- y)5p{y) (2) 

where the ideal gas contribution is proportional to In p and u is related to the direct corre- 
lation function. An apparently appealing feature of this model is that the ideal-gas contri- 
bution to the equation of motion satisfies 

v(pV^p)=TVV (3) 

There has been considerable effort to study the SDK system^. The most complete analysis 
is due to Kim and Kawasaki 16 . Analysis of the SDK model field theory uncovers multiple 
complications such as nonpolynomial nonlinear structure and multiplicative noise structure. 
This leads to an impractical nonlinear form of FDT relating response and correlation func- 
tions. In the end, one is not definitively able to answer the question: Is there an ENE 
transition at one-loop order for the SDK model? 

Our second order self-consistent microscopic theory suffers none of these problems and 
gives a definitive answer to the ENE transition question. An interesting point is how the 
theory handles convergence of integrals in perturbation theory. The theory naturally orga- 
nizes itself into a structure with self-dressed propagators. This leads to convergent integrals 
in perturbation theory. One finds that self-dressed propagators occur order to order in per- 
turbation theory. The dressing consists in multiplication of the physical correlation function 
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by factors of and the effective potential. It is crucial that one recognize that these 
self-dressed correlation functions themselves satisfy a FDT. 

It is extremely useful that the linear FDT holds in the theory This facilitates the use 
of a simple kinetic equation in treating the time evolution and separating out the static 
behavior. 

While not immediately apparent, the understanding of the role of one-particle irreducible 
vertex functions in this analysis is very important. The expansion parameter is not, as in 
conventional field theory, in terms of the vertices. There are nonzero nonlinear vertices even 
at zeroth order in the interaction. An example is the cubic 1PI vertex. The one-particle 
irreducible vertices are interesting even in the noninteracting limit. Starting at the four-point 
vertices one must deal with one-particle reducible contributions to the cumulant structure. 
Vertices are now generally frequency dependent. We reserve the discussion of the four-point 
vertices to the next paper in this series. 

The static theory can be developed in complete analogy with the dynamic theory. In 
particular one can work out the self-consistent static expansion in terms of the effective 
potential. One finds exactly the same results for the statics using dynamics and statics. 

The approach taken here is fully microscopic and allows one to calculate in detail. This 
is in contrast with the projection operator technique which allows one to develop a useful 
phenomenological description of a problem. There was a substantial effort 17 to show that a 
noninteracting Brownian gas can be described by a nontrivial but solvable cubic MSR field 
theory. Our work has been compatible with these results finding agreement for the set of 
Brownian-density cumulants. 



Let us begin by defining the dynamical system of interest. Consider a set of N particles 
with configurations specified by the coordinates Ri which satisfy the equations of motion 



II. SMOLUCHOWSKI DYNAMICS 



Ri = DFi + rji 



(4) 



where the particles experience force 







U(R), 



(5) 



dRi 
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with total potential 

u(R) = lj2v(Ri-Rj) (6) 

where V is a pair potential. There is a noise source rji for each coordinate which is taken to 
be Gaussian with variance 

(7*(*)tfc(0) = 2k B TD5(t - t')5 t] (7) 

where D is a diffusion coefficient. 

We showed in FTSPD that one can set up a useful self-consistent perturbation theory if 
we develop the theory in terms a small set of collective variables For this system, the 
density p is essential since it governs the static equilibrium behavior and, from the point 
of view developed here, is always included in $ = (p, . . .). The set $ must also include a 
response field B (described below) if we are to control and manipulate the interactions in 
the system. The set of collective variables treated [$ = (p, B, . . .)], is flexible and controlled 
by pairing each observable with a conjugate external field [H = (H p , Hg, . . .)] Here we 
specifically treat fluctuations in equilibrium and choose [$ = (p, B)]. We assume the system 
is in equilibrium initially [t = t , Ri(t ) = Rl ] and the initial distribution for a set of N 
particles is canonical: 

P [R°] = e-W R ^/Z (8) 

where U is the potential energy defined by Eq.(|6]) and j3 is the inverse temperature. The 
MSR action for the problem is given by 

N 

dti 

'to , i 



fOO N 

A R = dhJ2 IMtJkBTDRiitx) + iRiih) ■ (Riih) - DF^))} + Aj (9) 

J tn -i 



where the contribution to the action Aj is from the notorious Jacobian^. 

It was shown in FTSPD that this action can be written in the highly compact form 

A = A + A T (10) 

where A is the quadratic part of the action excluding the quadratic contribution to the 
initial probability distribution, 

N 



A = dh^RiDRi + iRi-Ri , (n) 
Jto i=l J 



where D = ksTD and where the interaction is given by 

A i = \Y.I ^2^(1)^(12)^(2). 
The Greek labels range over p and B and we introduce the interaction matrix 

<7 a „(12) = -/9F(12) [p a/ M(fi - to) - (f>*B v + B aPu , 
where we have introduced the useful notation 



Pa &api B a 5, 



aB 



and 



V(12) = V(x 1 -x 2 )5(t 1 -t 2 ) 



The conjugate field is given by 



N 

B(l) = D [(£iiVi + 0(0) V?)] 6( Xl - Riih)). 



i=l 



The canonical partition function can be written in the convenient form 



J N 



N 

n 

i=i 



P (RW)e- A °- A ' +H -* 



= Tr {N) e' Al+H ^ ', 



where we have introduced the average 



«=i 



V{R i )V(R i )d d R. 



■d p(0) 



P ( J R; o) )e- Ao C»(i?). 



In the grand canonical ensemble, 

oo N 

Z T [H] = J2 —\Tr [N) e$ dlH ^>^ e ^ <* 1 * 2 *( 1 )-^*( 2 ) > 



Af=0 



with the generator of cumulants given by 



W[H] = \nZ T [H}. 
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III. SELF-CONSISTENT DEVELOPMENT 



It was shown in FTSPD that the one-point average 

_6_ 

m 



Gi = = 4f W [H] (22) 



satisfies the identity 

d = Tr^e^ +Aw ^, (23) 
where i labels space, time and fields p or B, and where 

AW[H] = W[H + F] — W[H] (24) 



with 



F = J2^A (25) 



J 



and 



N 

^ = E# } - (26) 



a=l 



We have 



= *[(*i-i2< >(fi)] (27) 



and 



^ 0) (1) = D^fa)^ + 0(0)^]^ - flWfa)]. (28) 

These results were established in FTSPD using functional methods. 

The dependence of the theory on the interaction potential is controlled by the quan- 
tity AW [if] = W[H + F] — W[H]. We can expose the dependence on the potential by 
constructing the functional Taylor-series expansion 

AW [H] = N; F^Wm + N; ^F^Wm + . . . (29) 

I L %J L J 

and we can conveniently introduce the set of cumulants 

a ^ = -km--k wm (30) 
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to obtain 

AW[H] = £ FiGi + £ + E ^ W*GW + • • • ( 31 ) 

i ij ijk 

with given by Eq. (!25|) . Clearly, in this form we can take AW to be a functional of Gi. 
One can then use functional differentiation to express higher order cumulants in terms of 
the one- and two-point correlation functions Gi and Gij. 

Of particular interest is that we established in FTSPD a dynamic generalization of the 
static Ornstein-Zernike relation^. Taking the functional derivative of Eq. fl23j) . we have 

n - n 
13 6H, ' 

= Gij + E C ikGkj (32) 

k 

where 

Gil = TrMj^ AW (33) 

is a single-particle quantity and we have the memory function^, self-energy, or dynamic 
direct correlation function given by 

c l3 = fr^^^-AW. (34) 

Since AW can be treated as a functional of Gi we see at this stage that we have available a 
self-consistent theory. If we define the matrix inverses 

E ^ikGkj = Sij (35) 

k 

and 

E likGkj = Sij (36) 
k 

then the two-point vertex is given without approximation as 

Tij = jij + Kij, (37) 

where 

Kij = - E %kCkj (38) 

k 

is the collective contribution to the self-energy. 
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IV. FLUCTUATION-DISSIPATION THEOREMS 
A. Introduction 

The theory developed in FTSPD is very general and applicable to a wide variety of 
nonequilibrium problems. Here we look at fluctuations in equilibrium and see that one has 
a fluctuation-dissipation theorem (FDT) available. The existence of a FDT allows one to 
organize the theory very efficiently 

B. Time-reversal symmetry 

Let us focus on the time reversal transformation given by^ 



rRi(t) = Ri(-t) 



(39) 



rRi(t) = -Ri{-t)+ipFi(-t). 



(40) 



How does the action 



A 



/oo 
-oo 



■^Ri(t) - DFi(R) 



(41) 



change under the transformation? We have 



A' = tA 



/OO 
-oo 



[-Ri(-t) +iPF i (-t)}Dp- 1 [-R 4 (-t) +i/3Fi(-t)] 



+ i[-Ri{-t)-0Fi(-t)] 



^Ri(-t) - DF(tR) 



+ tA,, 



We need to look at the various terms. Consider first 



(42) 



/oo 
-oo 



ifiFii-t)] 



d_ 

Of 



Ri(-t) 
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/OO ft f—OO f) 

Next we collect the two force terms 

/OO 
dt2/3[(-i)i ? iH)^(-*) + ^(-*)^(-*)] = o 

Putting this together we have 



(43) 



A'=l dtJ2 

-oo 



[-Rii-WDp-^-Rii-t) + 2im(-t)} 



+ i[-Ri(-t) 



^Ri(-t) - DFi(rR) 



(44) 



where since Aj is a function of the density, tAj = Aj. Then combining the terms linear in 
F and letting t — > — t we have the invariance principle A' = A. 

C. Application of the invariance principle to averages of fields 

The first application of the invariance is to correlation functions. For the n-point density 
correlation function, we have, from its invariance under r, 



G P p... p (l, 2, . . . , n) = G pp ...p(l, 2, . . . , n) 



where t = (x£, —tf). 

Consider the two-point response function 



G pB (12) = (p(l)S(2)> 



(45) 



<p(l)B (2)> + ( P (l)Bj(2)) 



where the particle density is given by 



N 



i=l 



The response field is the sum of 



N 



B (2) = D^iRifo) -V X2 5[x 2 - Ri(t 2 )] 



i=i 



(46) 



(47) 



(48) 
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and 

Bj{2) = 6(0)DV 2 X2 p(2). (49) 

Applying r to B 

N 

t£? (2) = Dj^i-Rii-h) + iPFi(-t 2 )] ■ V X2 i5\x 2 - R t (-t 2 )] 



i=i 



= -B (2)-8DV X2 F(2). (50) 

We can then write 2 — (x 2 , — t 2 ) and 

r£(2) = -5(2) + 25,(2) - f3DV X2 F(2). (51) 

Then we have on applying the invariance principle 

G>(12) = -G pB (l, 2) + 20(O)DV 2 X2 G pp (l, 2) - f3DV X2 G pF (l, 2). (52) 

We now need an independent expression for the density-force density correlation function. 
As a first step construct an operator Bo (2) that satisfies 

roo N 

O B0 (2) dtY / R 3 (t)DR J (t) = B (2). (53) 

J-oo j=1 

It is easy to see that 

N iB 5 

Ob {2) = ^ fV* ■ jj—Ax, - R 3 (t 2 )}. (54) 
j= i 2 SRjih) 

We now use the identity 

Trf(R)6 Bo (2)e- A = (55) 

or 

(f(R)6 B() (2)A} = 0. (56) 
Letting the operator act on A, we have 

6 Bo (2)A = Bo(2) - ~P{2) + ^V X2 F(2) (57) 



m + Uw 2 ~ Dvl ) p(2) + ir v - F ( 2 )- 



Multiplying by the density and averaging, we have 
'p(l) 



B{2) + Udf 2 ~ Dvl ) p{2) + ^T V ^ F{2 ^ 



or 



G pB {l2) + 



0, 



DV^ G pp (12) + ^V, 2 G pF (12) = 



2\dt 2 ~ -**)-r^—> ' 2 
Eliminating G pF between Eqs. (l52|) and (I60j) we have 



G pB (12) = -G pB (l, 2) + 20(0)DV 2 n pp (l, 2) + 2G pfl (l, 2) 



- + G pp (l, 2) 



= G p b(1,2)-/3^G pp (1,2). 

Since G p # is retarded we can multiply by 9(ti — t 2 ) and obtain the crucial FDT 

G pB (12) = -t 2 )/3^G pp (l,2) 

ot 2 



^(ti-t 2 )/3^G pp (l,2). 



D. FDT: Fourier Transform 



The time Fourier transform of Eq. fl62|) is given by 



/oo >5 
Jte™ 4 pe(t)-G pp (q,t) 



o 



OO ,C) /-OO rjfi 



dt J-oo 2n 



/oo Aii roo 
— (-iu)G pp (q,u) / dte<«~ a * 
-oo Z7T JO 
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f°° du . ._. . _. % 
J-oo 2n u — ui + ir) 



So we have 



GpBM= rf m £f^L. (64) 

J-oo 2,71 uj — oj + m 



CO — CO + IT) 

Taking the imaginary part 



J-oo Z7T W — W + 



— ^ G pp(9> w )- ( 65 ) 



It is useful to check the FDT in the noninteracting limit. 



E. FDT and vertex functions 

Dyson's equation, 

is a matrix equation relating matrix elements for the cumulants to the matrix elements of 
the vertex functions. Using Gbb = and T pp = we have 

G p b = (67) 
Gb p = (68) 



and 



G PP = -^r BB -L. (69) 

J- Bp J- pB 



Since G p b = G* Bp and r Bp = r^ p , we have 

°pp = " (70) 

I 1 Bp] 

and since G pp is positive, we have 

Tbb < 0. (71) 
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We can now use the FDT in the following way. Taking the imaginary part of G pB , 

Im G pB = -Im ^ pB 

|t Bp | 

= ~^G„, (72) 
and using Eq.(l69l for G pp , we can cancel the positive denominators and find 

Im T pB = ^-^bb (73) 



or 



Im T Bp = ~~Tbb- (74) 



F. Kinetic equation and the FDT 

Starting with the Bp component of Dyson's equation in q, t space, we have 

J dtT Bp {q,t-t)G pp {q,t-t') + J dtr BB {q,t-t)G Bp {t-t') = Q. (75) 
Let us now write 

r a p(q,t) = 1 ^(q,t)-Z a p(q,t) (76) 

where 7^3(9, t) includes the zeroth and first order contributions in perturbation theory and 
the self-energy E starts at second order. Suppressing the wavenumber label, we can rewrite 
Eq. (!75l) in the form 

7B>)Gp P (M') +^Ut)G Bp (t,t') = tf(t,f), (77) 

where 

= f dsY, Bp (t - s)G pp {s - t') + f dsY, BB (t - s)G Bp {s - t') (78) 
J— 00 J— 00 

using the fact that S^ p (t — s) ~ 6(t — s) and G Bp (s — t') ~ 6(t' — s). We then use the 
fluctuation-dissipation theorems 

Z Bp (t-s) = 9(t-s)/3^Z BB (t-s) (79) 
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G Bp (s - f) = 6{t> - s)/3^-G pp (t' - s) (80) 



dt 

to obtain 

d 



m(t,t') = -f ds (3^ BB {t-s) G pp (s-t')- f dsE BB (t-s)p^-G pp (t' - s). (81) 

J — oo OS J—oc OS 

If we integrate the first integral by parts, we have 

tt(t, *') = -^ BB (0)G pp (t - t') + (3 dsH BB (t - s)^-G pp {s - t') 

ft' d 
-J dsT, BB (t- s)—G pp (t' - s) 

= -Z BB (0)G pp (t - t') + J* dsZ BB {t - s)^-G pp {t' - s) (82) 
where we assume t > t'. We then have the kinetic equation 

h { B l{t)G pp {t,t') = -p 2 Z BB (0)G pp (t - t') + /3 2 £ dsE BB (t - s)^G pp {t> - s). (83) 



We see that our dynamical problem is now in the form of a memory function equation 
and the dynamic part of the memory function is given by the self-energy T, BB . With some 
additional work, one can show the equal-time quantity 

/3 2 P o W<7, t = 0) = S-\q) - [1 + Po f3V(0)], (84) 

where V(q) is the Fourier transform of the potential. We will discuss these results in great 
detail elsewhere. 

G. Static Implications 

The integral form of the FDT tells us that in the small u limit 

/OO AHj 
— (-P)G pp (q,Cu) (85) 
-oo Z7T 

= -PS(q) (86) 
where S(q) is the static structure factor. One can also write this as 

S(9) = "OT (87) 

This result will be extremely useful to us. 
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V. IRREDUCIBLE VERTEX FUNCTIONS AND BROWNIAN GAS 

CUMULANTS 

A. Irreducible vertex functions 

We list here the fundamental definitions relating the cumulants in the theory to the 
reducible vertex functions we need. The two-point irreducible vertex is defined by the 
Dyson's equation 

FikGkj = (88) 
The three-point one-particle irreducible vertex is defined by 



which is equivalent to 

Gijk ~FTj~ Cri u Gj v Gk w T uvw (90) 

Ott k 

which can be rewritten using Eq. (IBB]) as 

^ijk j v Y ku>G U vw (91) 

The case of four-point cumulants and irreducible vertex functions is discussed elsewhere^. 

B. Noninteracting Brownian particle cumulants: wavenumber and time regime 

In FTSPD we derived a compact expression for the n-point cumulants for the fields 
<p = (p, B). In the time wavenumber regime, we have 

G { Za 2 ,...,a n (1, 2, . . . , n) = Po%i + <72 + • • • + q n )hh . . . b n e N - (92) 
where bj = b aj (qj,tj) with 

&„(<&,*,•) = !, (93) 



b B (q j ,t J ) = (3 ]T K jk 6(t k -t 3 ) } (94) 
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1 - 

N n = 7^H KalU-tjl (95) 

z i,i=l 

and 

#y = Dq % ■ qj . (96) 

This expression is manifestly translationally invariant in space and time. 
For our purposes, we need these cumulants Fourier transformed over time: 

G ( S,a 2 ,...,a n (h 2,...,n)=p 5(q 1 + q 2 + ... + q n ) 

x J dt x dt 2 . . . dt n b x b 2 . . . b n e^=^ kih e Nn , (97) 

where we have not found it confusing to use the same symbol for the Fourier transform. In 
this paper we assume the system is in an equilibrium state for all time. 

C. Time Fourier transform 

We discuss the method of taking these time Fourier transforms elsewhere^. Here, we 
simply list the needed results. We will need the zeroth order two-point cumulants 

Gg(12) = 2k iPq G 1 G\8{1 + 2), (98) 
Gg)(12) = -/3K lP oG^(l + 2), (99) 
G%{12) = -f3K lPo G 1 5(l + 2), (100) 



and 



where 



and 



Ggi(12) = 0, (101) 

Ki = Ku = Dql (102) 

G7 1 = -iui + K h (103) 

5(1 + 2) = (27v) d 5( qi + g 2 )27r<5(wi + w 2 ). (104) 
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D. Two-point vertices 

We also need the zeroth order two-point irreducible vertices defined by 

lS?GS?=V (105) 

We find 

/3 7 £>(12) = 0, (106) 
/3 7 g(12) = -^(1 + 2), (107) 

/3 7 g(12) = -^*(l + 2), (108) 

and 



^(12) = -— 5(1 + 2). (109) 
These are the key building blocks of the theory. 

E. Three-point vertices 

A key role in the second order theory is played by the zeroth order three-point irreducible 
vertices. These are defined in terms of the three-point cumulants 

T^(123) = -72 ) (l)7j?(2) 7 ffi(3)G tttw (123). (110) 

These can be found to be given by 

7£(123) = 7 £Ul23) = 0, (111) 



7 g>(123) = -^ K 12 G^'* + K 13 G^'* 6(1 + 2 + 3) = ^[l-iE 1 ]6(l + 2 + 3), (112) 



1 — 1,* i TS — 



1 



Po 
and 

(0) 1 ■- 



/3 7 ^ p (123) = -2-K 12 5(l + 2 + 3) (113) 
Po 



where 

5(1 + 2 + 3) = (2n) d 5( qi + q 2 + q 3 )27r5(co 1 + co 2 + cu 3 ) 

and 

Ei = u 2 K 12 + uj 3 K 13 . 
The other vertices can be constructed using symmetry. 

VI. FIRST ORDER RESULT FOR K 

Going to first order in the potential we have 

AW^=Y,F U G U . 

u 

We next need to compute 

which goes into Eq. (1341) giving the result 

U 

Putting this result into Eq.(l38l) and using Eq.(l36l) gives the very simple result 

K\f = a„. 
These results satisfy the FDT in a trivial way: 

K% = K§ = -0V(q) 

and 

K^ B = - A/ m K $) = 0. 

pui 
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VII. FULL FIRST ORDER SOLUTION 

In FTSPD we worked out the first order theory concentrating on the two-point cumulant. 
Here we work things out including the equation of state and one-point averages. This sets 
the stage for the second order calculation. 

At first order in perturbation theory^, we have the equation of state 

Gi = Gf ] + TrfaAW^ + ... (123) 
and for the two-point cumulant, the kinetic equation 

Gij = G?> + fr^AW^ + G\fa £u G uj . (124) 
First, we rewrite the first order static term 

= frfcAWW = fr&FuGu 

= fr^i<j>iat u G u = G[fa iu G u (125) 

and the two-point function 

= frtfiifafaaeuGu = Gf^a tu G u . (126) 
Thus, one has the coupled set of equations 

G% = Gf^ + G^ ai u G u (127) 

and 

Gij = G^ + G^l<7£ U G u + G$<jf u G U j (128) 
where in the field theory protocol 

a ij (k) = V(k)[ Pi B j +B iPj ]. (129) 
Let us look first at the term with the cubic cumulant 

g\f = G { °la eu G u . (130) 
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We can write 

r (0) _ r (0) r (P) n (P) (0) /i 01 \ 

^ijl ~ ^iu ^jv ^Iw luvw \ loL J 

and thus 

Mi) _ r (o) r (o) r (o) (o) r 

G{ u Gjv luvw (G w Cui"*) (132) 

where in the last step we used the equation of state, Eg. (I127[) . If we use the useful identity 

(o) G (o) = _ (o) (133) 

and realize that we can write 

Gi = ^-Gf ] (134) 
Po 

where p — (p), we have 

Q^^-G^G^i-^ip/po-l)} 

= 4 0) (P/Po-l). (135) 

The two-point equation becomes 

Gij = G^ + G\®\p/po — 1) + G^aeuGuj 

— + G^ciiuGuj. (136) 

Po 

Now p/ po is determined by the equation of state. If we assume the one-point average is of 
the form 

G u (l) = PuG^q^ixbM (137) 

then we have 

Pi[G p ( qi ) - G^(q 1 )}2n5(co 1 ) = Gf {q 1 ,co 1 ) ( j eu {q 1 )G u {q 1 ,uj 1 ) 
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= G ( $\qi,0)a eu (q 1 )p u G p (q 1 )2ir5(uj 1 ) 

= G^\q 1 ,0)a ep (q 1 )G p (q 1 )2n5(co 1 ) 

= G ( S(q u 0)V(q 1 )G p (q 1 )2ir5(u Jl ) 

= p i G%{q 1 ,0)V(q 1 )G p (q 1 )27r6(u 1 ) 

= p l (-p )/3V(q 1 )G p (q 1 )2n5(uj 1 ). (138) 
Canceling common factors gives 

G p (q) = G p °\q)-p o pV(q)G p (q) (139) 



or 



G ^ = TT^W G " ,(9) ' (140) 

In the homogeneous limit 

Turn next to the two-point correlations which satisfy Eq. (11361) . Taking the pB matrix 
element, one gets a closed equation 

G pB {\) = (p/p )G%(l) + G%(l)V( qi )G pB (l) (142) 

which can be solved to give 

G pB ( gi )= - " r '"Z:Z)7 : (1-13) 



(p/Po)gjS(l) 
1 - G%{l)V{ qi ) 



PpKl (144) 



-zwi + Ki[l + y(gi)] 

where K\ = Dq\ and V(qi) = po/3V(qi). One can then solve for the density- density correla- 



tion function 



G PP (l) = (p/po)G$Kl) + GfJl)V{ qi )G Bp {l) + Gf B {l)V{q x )G pp {l) (145) 
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or 



G »M = i M*Lm„ ^ (iP/Po)G ( $(l) + Gf p {l)V{ qi )G Bp {l)) . (146) 
1 _ G B {l)V{qi) 



Putting in the results for the zeroth order correlations and 



leads to the final result 



GgJ(l) = (147) 



u( + K{[l + V(qi)} 2 



Turning to the statics, we can use the result found earlier 



/3S(q) = G pB (q,0) = (149) 

K[l + Vfa)] 



or 



s ^ - irk)' (150) 

Comparing Eq. (ll50j) with the static Ornstein-Zernike^ relation we can identify the effective 
interaction 

V EFF (q) = -p- X c D {q) (151) 

where Cd(q) is the physical direct correlation function which is assumed to be known by other 
means. We can, for example, assume that cjj(q) is given in the Percus-Yevick approximation 
for hard spheres 6 . 

VIII. SELF-CONSISTENT COLLECTIVE SELF-ENERGY AT SECOND ORDER 

The evaluation of the collective part of the self-energy, K, at second order begins with 
determining fj, defined by 

fi = 4r^W (152) 



SG j 



at second order. We have 



AW® = l -F u F v G uv (153) 
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and 



(2) 6 AW ( 2 > _ 1 8G UV 

f > -^g~-2 FuFv ig;- (154) 



We then have the standard functional manipulations, 



Sir* " *J / y ~ ™ r/"< " uv 

= —^2GiuGj V T^ k (155) 

Ml) 

where we introduce the three-point vertex 

T ijk = jTT^ij- (156) 
6G k 

Putting this back into Eq. fl 154ft gives 



p ( 2 ) SAW® I 
' j SG, 2 



fj — — — ~~ —-^FiF n G£ U G nv T UV j (157) 



This in turn goes into Eq. (j34j) and 



r^ 1 — ——Trrh F F C C V 



n Q its@ ln@ smG nuG m v^- uvj 



where 

£U = Tr <f) i( f>i<f) s e AW (158) 

is a three-point self-correlation. Things can be written more symmetrically if we introduce 
the three-point self-vertex 7^ via 

(159) 

Then, the collective part of the two-point vertex is given at second order by 



lit' s' Gil' Gss'O'iuO'sv G ur G vq T rq j 
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-'y g g r 



where we have the self-dressed propagators 

Guv Quv&vs 

G sr . (160) 

We see already the structure for making vertex corrections. At lowest order in the interaction 
we have 

rgi = 7$ (i6i) 

and we have the nontrivial approximation for the second order contribution to the collective 
part of the two point vertex (self-energy) given by 

k (2) - ~->y {0) r r -v (0) (it\9) 

ly ij 2 "<«» '"' "9 Ikrq- l xu ^V 

Clearly K can be constructed to be symmetric and therefore G can be constructed to be 
symmetric: 

) • (163) 

Equation (11621) is the key result to be analyzed carefully. We have a collective contribution 
to the two-point vertex which has the following desired properties: 

(i) It is quadratic functional of the physical density- density correlation function. 

(ii) We will show that it satisfies the FDT. 

(iii) In the short-time limit, it gives the correct static contribution to the static structure 
factor at second order in the effective potential. 

(iv) It determines the phase diagram for ENE transitions. 

In addition to the collective contribution to the two-point vertex, there is a single-particle 
contribution. This part of the self-energies is very interesting and will be fully treated in a 
companion paper. 



IX. SELF-DRESSING PROPAGATORS 



A. Second order self-energies 



It appears not to be a coincidence that the matrix propagator 



o ^ 1 @uvG V j ~\~ Gi U u uv G V j (164) 

uv 
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appears in the one-loop expression for the collective part of the second order self-energy. 
In treating the single-particle contribution to the second order self-energy we find another 
quantity 

Gij ^ ] G^ u ^uvGr V p<JpqG gj , (165) 
uvpq 

which appears in the theory. These complications turn out to be welcome since G and G 
can be treated as effective matrix propagators which satisfy the FDT themselves. They also 
approach zero faster than G as q and cu go to infinity, thus ensuring convergence of integrals 
in perturbation theory. We expect additional similar quantities to appear at higher order in 
perturbation theory. 

B. G- frequency regime 

Consider the effective matrix propagator in Fourier space[(l = (qi,Ui)] given by 

G^(l) = \ E [^(1)^(1)6^(1) + G Q ,(1)^(1)GS(1)] . (166) 

We assume here that we are working in the field-theory protocol where 

<V(fc) = V(k)(pfj,B v + p u B p ). (167) 

Taking components of Eq. (11661) we see that the symmetrization does not influence the re- 
sponse contributions: 

G p b{1) = W [Gf»{l)a, u (l)G vB (l) + 0^(1)^(1)6^(1) 

= 5EK5(i)^(i)G pB (i) + G w (i), w (i)G; o ](i) 

= \ \G p %{\)V(\)G pB {\) + G pB {l)V{l)G p %{\) 

= G%{l)V{l)G pB {l). (168) 

G p b and G pB have the same analytic structure; they are analytic in the upper half-plane. 
Next, 

G BB = 0, (169) 
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while for the density-density component, 

= \ E [^(1)^(1)^(1) + G w (l)<v(l)Gg(l)] 

fJ,V 

= \v{l) [G%{1)G PP {1) + Gf p {l)G Bp {l) + G pB {l)Gf p {l) + G pp {l)Gf p {l) 

= \v(l) [Gg)(l) (G fl ,(l) + G pB (l)) + G pp (l) + 

= Gf p {l)V{l)ReG pB {\) + G„(l)V(l)iteG<g(l) (170) 

which is real but not necessarily positive. It is crucial to realize that G itself satisfies the 
FDT if G and G<® satisfy the FDT. 
The proof is as follows. Starting with 

GpaO) = Gf B {l)V{l)G pB {\) (171) 

and assuming V(l) is real, we take the imaginary part, 

Im G pB (l) = Im G%(l)V(l)Re G pB (l) + Re G p %(l)V(l)Im G pB (l). (172) 

Multiplying by — 2/(/3o>i), we have 

- ImG pB {\) = (im G%(l)V(l)Re G pB {\) + Re G%(l)V(l)Im G pB (l)) 

= Gf p {l)V{l)Re G pB {\) + Re G p °j(l)V \l)G pp (l) 

= G pp (l). (173) 



Thus G satisfies the FDT 



and we can write 



G p » = -^ImG pB (uj) (174) 

pui 



G pB {u) = I f- (1?5) 



and 



27T w — w + ii] 



r du PuG pp (u) ^ 

J 2lT id — Ld — in 



2tt uj — to — irj 
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C. G-time domain 



In the time domain we have the convolution 

G PP {t) = -V(k) / ds (G Bp (t - s) + G pB (t - s)) G p ° p \s) 

+ l v{k) SZ dsGpp{t - s) ( G ^ (s) + G ^ (s) ) • 

This together with the fluctuation-dissipation theorem leads to the short-time results 

G pp (0) = -(3V(k) Po S(k) 

G pp (0) = 

d pp (0) = pV(k)G pp (0)G p ° p \o) = 0. (177) 

X. SECOND ORDER COLLECTIVE SELF-ENERGY AND THE FDT 

We have the crucial one-loop contribution to the self-energy 

T% c \\2) = -i|d3d4d5d67^(134)G^(35)G^(46) 7 £ f (256). 
In terms of Fourier transforms in space and time we find 

rgf(Pi,P2) = i^) d+1 s(pi +P2W% C \ P i) (178) 

and 

ifA-Pi) = "/ dpsdp^>(13A)5( Pl +p 3 +p 4 ) 

x ^(3)^(4)7^(1,3,4) (179) 
where we have introduced the notation 



G is a correlation function that satisfies the fluctuation-dissipation theorem in the form 

dui (3uG pp (cu) 
2n uj — u + ii] 

and 



rdopoG^^ 
J 2n ui — uj — irj 
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A. T Bp 

We can break Eq. (jl79j) up into components and associate a set of one- loop contributions 
which differ by different vertices and propagators. A number of self-energy contributions 
vanish due to causality and for T Bp we obtain three contributions. (To simplify the notation 
we suppress the superscript on the three-point vertices.) Thus 

rg p C,1) (-Pi) = -~/ ^3^47^(134)5(^+^+^)^(3)^(4)7^(134), (183) 
rgf 2 \-Vi) = -\j ^3^47^(134)5(^+^+^)^(3)^(4)7^(134), (184) 



and 



rgf Z \-Vi) = -~/ dp 3 dprf Bpp {13A)8{p x +p 3 +p 4 )^(4)G pp (3) 7pAB (134). (185) 

We want to write these expressions in terms of G pp and we must exhibit care in treating the 
frequency integrals. Noting that /3 2 7 p bb(134) = — 2_K" 34 , we have 

0rgf 1] (-pi) = -\\ dp 3 dp,(3 lBpp (13A)5(Pi + p 3 + P4 )(3- l G pB (3)(3- 1 G pB (A)(-2K u ) 
= I dp 3 dp i (3 lBpp (l3A)5(p 1 +p 3 +p i )r 1 G pB (3)(K u ) [ ^ ^ G ^\ 

J J An U4 — + 17] 

= J dp 3 dp^K u J ( ^-u) A G pp (u 4 :)l3'y Bpp (13A)5(q 1 + k 3 + fe 4 )5(wi + w 3 + w 4 ) 



CO4 — + 277 



J dk 3 dkiK 34 5{qi + k 3 + fc 4 ) y ^^ 4 <5 w (w 4 



2tt y " ^ w -w 1 -w 3 -W4 + ^ 

dk 3 dk 4 K 3i 5{qi + k 3 + k 4 ) J -^-u)iG pp {u A ) 



31 



2tt 



x / ^/37b pp (134)/3- 1 G p b(3) 1 _ -. (186) 

Now G p b(3) is analytic in the upper half plane. Closing the contour in the upper half plane, 
£rgf ''(-pi) = (-1) / dhdhK u 5( qi + h + h) I ^^ 4 G pp (co,) {2m) 

x r'GM-ui - W4)/37b pp (1, "I - 4,4) 
= y dk 3 dk A K 3A 5(qi + k 3 + k 4 )(-iK 34 ) J ^ ^ ^^4^(0)4) 

x ^3^3) ^ (!,_!_ 4, 4) 

—Ui —U 3 — U 4 ) + 17] " 

= J dk 3 dk 4 K 34 S(q 1 + k 3 + k 4 ) ^ ^ ^ ^^/3G pp (u 4 )u 4 u 3 
x ^(^3)^34)^/37^(1,-1-4,4) 

= OVSll (187) 

where we introduce the notation 

6[J] = J dk 3 dk 4 S( qi + k 3 + h)J^J ^G pp (cu 3 )G pp (u 4 )J (188) 

and 

J<g = ^(^34)^-/37^(1, -1 - 4,4). (189) 



Next we notice that V%f' 2 \- Pl ) = rgf ,3) (-Pi), 



so 



/3rg P +3) (-Pi) = -4 / rfP3«7 Bw (134)5( Pl +p 3 + p 4 ) 

x /3- 1 G pB (4)G pp (3)/3 7 p^(134) 
= -J dk 3 dk 4 S( qi + k 3 + k 4 ) J dco 3 ^* Bpp (ll -1 - 3)G PP (3) 
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< / dS> 4 - . /37pp^(13, -1 - 3) 

J — Wi — OJa — UJd + IT) 



-OJl — U 3 — W4 + 17] 

where 

jg p +3) = w 4j R-/37b pp (13, -1 - 3)/3 7p pb(13, -1 - 3). (190) 
The total is thus given by 

J Bp = RJ^& z Kufo* B Jl, -1 - 4, 4) + oo 4 (3 lBpp (13, -1 - 3)f3 lppB (13, -1 - 3)]. 
We next need the vertices 

/3 7 b pp (1, -1 - 4, 4) = 1 + ^ - iQX 13 , (191) 



and 



/37b pp (13, -1 - 3) = 1 + i£x - iQK 14 , (192) 



/3 7ppB (13, -1 - 3) = 1 - *£ 4 , (193) 



where 



Q = ui + u 3 + W4 (194) 



and 



£ 4 = wii^u + u 3 K u . (195) 

We then have 

j Bp = R_[iu) 4 ui 3 K u [l + iE l - iflK 13 ] + w 4 (l + i-Ei - ifttfi 4 )(l - i^)]- (196) 

If we group together the terms proportional to fl, set R^Q = 1, and notice that inside 
the integrations each term vanishes because of odd frequency integrals over G pp , we then 
have 

J Bp = R-[ito 4 co 3 K 34 [l + iEi] + w 4 (l + i#i)(l - iE A )\ 
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= R- [1 + iEi] [iu 4 u 3 K 34 + 004(1 — iE A )] 
— R-W + iEilliuiUsKu + u> 4 (l - i(uiK u + £3X34)] 

= i2_[l + iSi][a>4(l-iwi^i4]- (197) 
Inside the integrals we are free to symmetrize with respect to 3 -H- 4: 

J Bp = + iE 1 }[^[u 4 + w 3 ] - Jy£i] 

= ii? [1 + -ui- iuxE^R- 



.[l + iEl ]-'^[l + iE 1 } 2 R_ 



and 



pT Bp (- Pl ) = + iSx] - y[l + i£i] 2 i2-]. (198) 

Since 0\\iEi] = we have 

/?rgf (-pi) = 0[i - ^[1 +^] 2 i?-] (199) 



or 



/3rgf (ft) = o[\ + ^iL[l + i£ 1 ] 2 ], (200) 

where 

R_ = - — (201) 

cos + 104 — u>i — i?7 

and 

0[J] = J dk 3 dk 4 5(- gi + k 3 + h)J^J ^G pp (Cu 3 )G pp (Cu±)J. (202) 
We need the result 

Im PT%f\ Pl ) = ^6[R.[1 + iE,} 2 }. (203) 
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B. T BB 

We turn next to the second order self-energy with two £Ts. After using causality, there 
are five nonzero graphs. The first is the simplest and given by 

1 r 



/3 2 rgf' 0) (-pi) = --J dp 3 «7^(134)5(p! +p 3 + P 4)G pp (3)G pp (4)(3 lBpp (lM), 
while there are four similar contributions given by 
/? 2 rgf' 1} (-Pi) = ~\J dp 3 dp 4 (3 lBpp (13A)5( Pl + ps +p,)r 1 G pB (3)G pp (4)p 2 lBBp (l34), 

P 2 r { BB 2 \-pi) = -\j dp 3 dp 4 (3 lB J13A)5( Pl + p 3 + p 4 )G pp {3)(5' 1 G pB {A)(5 2 lBpB {\3A), 

/3 2 rgf ' 3) (-Pi) = ~\J dp 3 dp^ 2 lBpB {l3A)5{ Pl +P3 +Pa)G pp (3)(3- 1 G Bp {A)(3 1Bpp {IM), 
and 

2 rgj?' 4) (-pi) = / dp 3 dp^ 2 lBBp {l3A)5{ Pl +p 3 +p,)(3- l G Bp (3)G pp (A)(3 lBpp (13A). 
Then, the simplest graph can be rewritten as 

/3 2 rgf' 0) (-pi) = a[jffl] (204) 

where 

4°l = -^(fi)| 7 Bpp(134)| 2 (205) 

= - R + )[l + E\\. (206) 

Turning to the other four contributions, they can be grouped together in a way that shows 
their sum is real 

^ 2 rg^' 1_4) (-pi) = J dpzdpiSipx + p 3 + Pi)G pp {3)K 1A 
x \ lBpp {\3A)G pB {A) + 7 Bpp(134)G B p(4)]. 



This leads to 



P 2 r^\- Pl ) = d[j^ 



BB 
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where 



,(1-4) 
J BB 



2UJ4K14 J duj^5{uji + co 3 + U4) 



[1 + iE,] [l-iE 1 ] 



CO4 — LO4 + if] 004 — 004 — irj 



= 20J4K14 



[1 + i#i] [1 - i^i] 



-Wi — W3 — W4 + 17] —U>i — LU 3 — OO4 — 17] 



= -2u 4 K u ^R_[l + iE x ] +R+[1- 



and 



(207) 



The terms generated from ^ proportional to give zero contribution after integration over 
u 3 and Co 4, so 



Combining the two contributions gives 

Jbb = i\{R-~ R+)[l + El] = - J E?i| J R_[l + iSi]+ J R+[l-i J E?i] 



-[l+^Ei] 2 ^--[l-^Ei] 2 J R + 



-7mi2_[l +iE 1 ] 



Finally we have the result 



^bb\pi) = -ImO[R-[l + zE,} 



and we see that this set of self-energies satisfies the FDT as 



(208) 



(209) 
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XI. STATICS AND AN EFFECTIVE POTENTIAL APPROACH 

Let us summarize our basic results up to now. The two-point vertex is given in the form 

r^T^ + ^ + rg^TT + rr (210) 

where 7^ is the noninteracting gas result, 7^ and are the first order single-particle 

and collective contributions (respectively) to the self-energy, and 'Ji]'^ and are the 

single-particle and collective contributions at second order in perturbation theory. Most of 
the focus in this paper is on the collective part of the two-point vertex at second order. We 
will be able to show that the collective part determines the static structure to this order and 
determines the ENE phase diagram. The single-particle contribution does not enter into 
either determination. However, we anticipate that it plays a crucial role in the slow kinetics 
near the ENE transition. It will be treated in detail in a companion paper. 

Let us now determine the static structure in the second order approximation. This is very 
conveniently done in the case where there is a linear fluctuation-dissipation theorem. The 
key result is that the zero frequency limit of the response function gives the static structure 
factor 

G pB (q, 0) = - / %G pp (q, u) = -0S(q) (211) 
in the ergodic phase. Then, using the Dyson's equation we have 

-l/(/3p ) -V{q) + T Bp (q,0) 



or 



1 + pol3V{q) - /3p T Bp (q, 0) 



The needed self-energy is given by Eq. ( 12001) as 

f3T Bp (q,0) = f3O 



= A-{ dUA{2-K)%- qi + k 3 + k±)G pp (3)G pp {4) 

h I §7§? m > + + *•> I ^ (3) / ^ (4) - (214) 
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We showed previously that 

/ ^G pp (3) = G pp (k 3 ,t = 0) = -p f3V(k 3 )S(k 3 ), (215) 

J 271 

therefore, 

PT Bp (q,0) = ^J 0- d ^ d (2K) d 5{ qi + h + k 4 )p pV(h)S(h)p pV(h)S(h). 

Then, after introducing dimensionless wavenumber k' 3 = k 3 a (where a is the hard-sphere 
diameter), S(q') = S(q)/p Q , V(q') = p f3V(q), and 

S(q') = , , (216) 

V ; 1 + V{q') - M{q') V ; 

where 

M(g') = p pr Bp ( q , o) = ^ / + ^ + fcyvc^^v^)^, (217) 

where 1] = 7rpa 3 /Q is the packing fraction for hard spheres. For simplicity we drop the primes 
on the dimensionless wavenumbers. 

So far the analysis has been rather general. Let us apply these results to the case of 
three-dimensional hard spheres. We immediately have a problem if we view our perturbation 
theory expansion in the conventional way since the Fourier transform of a hard-core potential 
is not well defined. We can however take a different approach. Instead of taking V(q) as 
given and computing S(q) , we take S(q) as given and, at a given order, determine the 
effective potential by inverting the equivalent of Eq. (I216[) . Thus, at first order 

= ttW (218) 

and the effective potential is essentially the direct correlation function given by 

V(q) = J--l. (219) 

S{q) 

We assume that the structure factor for a hard-sphere system is given by the solution^ 
to the Percus-Yevick equation. In Fig. 1 we plot the structure factor for different packing 
fractions. 

The numerical solution of Eq. (I216j) for V(q) is obtained as follows. First, rewrite Eq. (12 16ft 
in the form 

V = S- 1 - 1 + M(V). (220) 
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Next, define 



V = S- 1 - 1 (221) 



and rewrite 



V = V + Q(V), (222) 



where 



Q(y) = M [V + Q{V)\ = 7(g). (223) 

We seek an iterated solution. To this end, we write 

<xQe+i + (1 - a)Qe = I(Q e ) (224) 

with the number a chosen appropriately. This equation can be rewritten as a functional 
recursion relation, 

Qt+i=Qt + -[I(Qe)-Qe]- (225) 
a 

At a fixed point we have a solution to the original problem. By trial and error we find that 
a reasonable choice for a is a = 100. Then beginning with Q = we monitor 

A, = ^ q (V e+1 -V e )* (226) 

Over the entire range of 77 studied, we find that is driven to small values (10~ 6 -10~ 8 ) 
with a minimum controlled by the range of wavenumber included in the calculation. One 
obtains good accuracy if one chooses q max = 80. 

In Fig. 2 we give the results for the effective potential for a sequence of densities. Also 
shown is V . A key quantity is VS. It is plotted in Fig. 3 where we also plot V S. 
Clearly the two quantities are close and do not change dramatically in going from first to 
second order. The biggest shift is for small wavenumbers. That M(q) has a maximum for 
small q is understood by assuming that V{k)S{k) is sharply peaked for k = k . Therefore 
V{k)S(k)V{k + q)S{k + q) is small except when k \q + k\ ~ k which is true for small q. 
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XII. ENE TRANSITION 

As we increase the density, (or lower the temperature), our model system slows down. 
It is possible that there is a density above which the density- density correlation does not 
decay to zero for long times. In the nonergodic phase 

limG pp (q,t) = F(q)>0. 

In terms of the time Fourier transform, 

G pp (q,u) = F(q)27r5(Lu) + g pp (q,u), (227) 

where g pp is the regular contribution. In the nonergodic phase, assuming the FDT still holds, 
one has for the response function 

/oo 
G pB ( q ,t) 
-oo 

= fi J^J(t)^G pp ( q ,t) 
= P J o °° ^G pp (q, t) = p[F(q) - S(q)} (228) 

while 

/oo 
G Bp {q,t) 
-oo 

= P J™J(-t)(-l)^G pp (q,t) 

= P f_J-l)jG pp (q,t) = -/3S(q). (229) 

There are several equivalent ways of determining the nonergodicity parameter F as a 
function of density and temperature. One approach developed in the time-domain results 
from pursuing the memory-function equation, Eq. (I83|) . This will be explored elsewhere. 
Here we give a more direct complementary analysis in frequency space. 

Using one of the Dyson equations we have 

G pp (q, u) = -G pB (q, uj)T BB (q, u))G Bp (q, u). (230) 
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For an ENE transition we associate a zero- frequency delta-function in V bb 

r BB (q,u) = —2ii5(u)T(q) + regular part. (231) 
Then, using Eqs.fl22B and f[23Tj) in Eq. f[230l . we find 

F(q) = p 2 T(q)[S(q)-F(q)]S(q). (232) 



Equation (12321) can be rewritten as 

F(q) 



S(q) - F(q) 



S(q)(3 2 T(q). 



We need to extract the lowest order approximation for T. We have from our previous work 
that 

/3 2 rgkl)= ImO[R-(l + IE,) 2 } 

x ImO[R-(l + iEi) 2 }. 
In the nonergodic phase we have the ^-function contribution 

G pp (k 3 , w 3 ) = -p f3V(h)F(h)27c5(u 3 ) (233) 
and we find immediately that we can set E\ — 0, and 

ImR- = 7r5(wi + oj 3 + u)i) = 7r5(u>i). (234) 
Comparing with Eq. (I231j) . 

(3 2 T B (q) = i / ^ Po /3V(fc)F(% /3^(g + A ; )F(g + A;). 
It is useful to write 

F(q) = S(q)f(q) (235) 
and it is /(g) that is conventionally called the ergodicity parameter. 
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Using the same set of dimensionless variables as in the static calculations, the nonergod- 
icity equation can be written in the conventional form 

where 

W{q) = Tk §iq) I + ^ + ^)S{q + k). (237) 

Before solving Eq. (I236j) numerically, it is useful to obtain an approximate analytical 
solution. Notice that the combination V(k)S(k) (see Fig. 3) is sharply peaked at the 
structure factor maximum. Then, to a reasonable approximation, we can replace f(k) with 
its value at the maximum such that 

w(q) = f 2 (q )S(q)M(q) (238) 

with M(q) the same quantity that appears in the static calculation. Putting this result in 
the nonergodicity equation we obtain 

f(q) 
i - /(?) 

where 



f 2 (qo)w(q), (239) 



w(q) = S(q)M(q) (240) 

is known from our static structural calculations. We now have a closed algebraic equation 
for /(go) if we set q = go m Eq. (12391) with f = /(go) and wq = w(qo): 



f0 r2 



f£w Q . (241) 



l-/o 

Discarding the ergodic solution / = 0, we have a quadratic equation to solve given by 

wofo - w fo + 1 = (242) 

with the solution 

1 



/o = 2(l + Vl-4M,j. ( 243 ) 

Note that there is no transition for wq < 4. For wq > 4, we obtain the full wave number 
dependence by putting / back into 
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It is easy to compute the Wq numerically using the Percus Yevick (PY) structure factor. 
In Fig. 4 we plot wq versus rj. We find that the ENE transition density, when wq = 4, is 
rf = 0.53. The associated f(q) is shown in Fig. 5 for several densities in the nonergodic 
phase. 

Returning to the full problem, Eq. ( 12361) can be rewritten in the form 

The solution of this equation is available via direct iteration. If we use an initial trial value 
°f foil) = 1/2 in Eq. (j237p . we generate an initial trial value w = 0.25S(q)M(q) to obtain 

fi{q)=w /(l + w ). (246) 



We continue iterating via 



fi+i(q) = w t /{l + w t ) (247) 



and find that 



lim (£+!-/*) = (). (248) 



i— >oo 



For pure hard-sphere systems in three dimensions we find an ergodic-nonergodic transition at 
rj* = 0.76. Clearly this density is physically unattainable. It is interesting to solve Eq. 



using the first order result V — S^ 1 — 1 rather than the second order result. Despite the 
fact that there is small change in the effective potential in going from first to second order 
one finds a substantial change in the critical density from 0.76 to 0.60. Thus, we see that rj* 
is a sensitive quantity. 

In Fig. 6 we plot f(q) for hard spheres for a set of rj > rj*. We find a lot of structure in 
/(<?). Comparing the approximate f(q) given in Fig. 5 with the full numerical solution we 
see good agreement despite very different transition densities. 



XIII. CONCLUSIONS 



We established here that the theory of time-dependent fluctuations in equilibrium can 
be organized around the fluctuation-dissipation theorem. Using the FDT one can relate 
the two-point response function to the two-point correlation function. In turn this gives 
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a linear-relation between self-energy components. We show more specifically the nontrivial 
result that the collective parts of the self-energies at second order satisfy a FDT. We find the 
second order self-energies as functional of the physical density-density correlation functions. 
The role of the three-point vertex functions are crucial in this analysis. 

Using the FDT the static structure is separated from the dynamical information and 
gives one a self-consistent expression for the static structure factor in terms of the potential. 
This expression can be recovered from a strictly static analysis. We show how the problem 
can be turned around and posed as in MCT. The structure factor is assumed to be given 
and we solve for the effective potential V that produces the known structure factor. At first 
order, the effective potential is proportional to the static direct correlation function. At 
second order, V is similar to the first order result except at small wavenumbers. 

We also show that the theory is compatible with an ENE transition at second order. 
The critical density is sensitive to the details of the calculation. Using the Percus-Yevick 
approximation for hard spheres we find an ENE transition at the physically inaccessible 
density rj* = 0.76. 

In a future work we will explore the role of the single-particle dynamics which occur in 
this theory. The treatment of these degrees of freedom gives one information about the 
equation of state, (p) = p(p ,T). The single-particle degrees of freedom contribute a second 
order contribution to the self-energies which depends on the full density- density correlation 
function. 

Our focus here, because of its simplicity, has been on Smoluchowski dynamics, but as will 
be discussed elsewhere, the method developed here can be applied to Newtonian dynamics 
and Fokker-Planck dynamics as well as a broader class of models. 

The approach presented here allows for a systematic method for analyzing corrections to 
this second order result including determination of higher order correlation functions. We 
plan to analyze the third order contribution to the collective self-energy soon. 

We guess that the theory presented here can be organized to give a theory of freezing. 
One can then do simultaneous free energy comparisons between the nonergodic state and 
the crystalline solid state. 



44 



Acknowledgments 

I thank the Physics Department and the Joint Theory Institute at the University of 
Chicago for support. A special thanks is due to M. Zannetti for showing that G p b satisfies a 
simple linear fluctuation-dissipation theorem and for other comments on this work. I thank 
David McCowan for help with the figures and comments and I thank Paul Sprydis and S. 
Das for their comments on this work. 



1 G. F. Mazenko, Phys. Rev. E 81, 061102 (2010). (Referred throughout as FTSPD.) 

2 A. Andreanov, G. Biroli and A. Lefevre, "Dynamical Field Theories for glass forming liquids," 
J. Stat. Mech.: Theory Exp. PO7008 (2006). 

3 S. Das, Rev. Mod. Phys. 76, 785 (2004). 

4 A. Crisanti, Nuc. Phys. B796, 425 (2008). 

5 W. Goetze in Liquids, Freezing and Glass Transition, edited by Hansen J. P., Levesque 
D. and Zinn-Justin J. (North-Holland, Amsterdam) 1991; and S. Das, Rev. Mod. Phys. 76, 785 
(2004). 

6 There is an exact solution to the Percus-Yevick approximation for hard spheres. See the discus- 
sion in Chapter 4 in J.-P. Hansen and I. R. McDonald, Theory of Simple Liquids, Third 
Ed., (Academic, New York 2006); N. W. Ashcroft and J. Lekner, Phys. Rev. 145, 83 (1966). 

7 A. Crisanti (private communication). 

8 The F\2 model and stretching are discussed in detail by W. Goetze, in Liquids, Freezing and 
Glass Transition, edited by Hansen J. P., Levesque D. and Zinn-Justin J. (North-Holland, 
Amsterdam) 1991. 

9 P. N. Pusey, J. Phys. Condens. Matter 20, 434202 (2008). 

The classical operator formalism for Newtonian dynamics was introduced by J. von Neumann, 
Proc. Natl. Acad. Sci. U.S.A. 18, 70 (1932); H. Koopman and J. von Neumann, ibid 18, 255 
(1932). 

1 The history of Smoluchowski dynamics is discussed in FTSPD. 

2 A brief discussion of the storied history of Brownian motion is given in Refs. 2 and 20 in FTSPD. 
We will refer to the noninteracting limit of Smoluchowski dynamics as the Brownian gas. 



45 



13 D. S. Dean, J. Phys. A29, L613, (1996). 

14 K. Kawasaki and S. Miyazima, Z. Phys. B: Cond. Matter 103, 423 (1997). 

15 MSR actions are discussed in Appendix A in FTSPD. 

16 B. Kim and K. Kawasaki, "A FDR-Preserving Field Theory for Interacting Brownian Parti- 
cles: One-Loop Theory and MCT," AIP Conf. Proc. 982 (Complex Systems: 5th International 
Workshop on Complex Systems," Sendai (Japan), Sept. 22-28, 2007.), 223 (2008). 

17 R. Velenich, C. Chamon, L. Cugliandolo, and D. Kreimer, J. Phys. A41, 235002 (2008). 

18 See the discussion in FTSPD. 

19 The static Ornstein-Zernike relation that connects the radial distribution function to the direct 
correlation function. 

20 In different contexts the kinetic kernels have different names. In the field theory context the 
kernels are typically called self-energies, in the kinetic theory case, where the analysis is in terms 
of retarded quantities, the kernel is called a memory function, and in the general case it can be 
called a dynamic direct correlation function. 

21 Note that b\ depends on the n U while the operator depends only on t\. 

22 G. F. Mazenko, D. D. McCowan and P. Spyridis, in preparation. 

23 An advantage of working in equilibrium is that we can take the initial time to minus infinity. 
We called this the field-theory protocol in FTSPD. 



46 




o 1 — 1 1 1 1 

5 10 15 20 25 30 

q 



(a) 



17.5 
15 
12.5 




5 10 15 20 25 30 

q 



(b) 







. 1 







5 10 15 20 25 30 

q 



(c) 

FIG. 1: The static structure factor as determined using the Percus-Yevick Hard Sphere Approx- 
imation for (a) rj = 0.1, 0.2, 0.3, and 0.4, (b) r\ = 0.5, 0.6 and 0.7, and (c) rj = 0.8. Larger first 
structure peak height corresponds to larger rj within each plot. 
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FIG. 2: The zeroth order potential Vo{q) and the numerical solution for the full potential V(q) 
from (a) 7] = 0.1 to (h) rj = 0.8 in increments of 0.1. 
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FIG. 3: The potential times the static structure factor for both Vo(q) and V(q) from (a) r] = 0.1 
to (h) r\ = 0.8 in increments of 0.1. 
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FIG. 4: Values of wo versus 77. Note that the line crosses wo = 4 (and therefore enters the 
nonergodic regime) near i] = 0.525. (The solid line is to guide the eye only.) 
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FIG. 5: Using the analytic approximation for f(q) given by Eq. (|244p . we plot the approximate 
form of the nonergodicity parameter for (a) rj = 0.53 (just above the analytic transition), (b) 
i] = 0.55, (c) r) = 0.60, (d) 7] = 0.65, (e) r] = 0.70 and (f) rj = 0.75. 
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FIG. 6: Numerically, we find a transition at 77 = 0.763. We plot the numerical solution and the 
analytic solution for (a) rj = 0.763 and (b) r? = 0.8. The more ragged and slightly lower line in 
each plot is the analytic approximation. 



